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We study the motion of a grain boundary that separates two sets of mutually perpendicular rolls 
in Rayleigh-Benard convection above onset. The problem is treated either analytically from the 
, corresponding amplitude equations, or numerically by solving the Swift-Hohenberg equation. We 

find that if the rolls are curved by a slow transversal modulation, a net translation of the boundary 
follows. We show analytically that although this motion is a nonlinear effect, it occurs in a time 
scale much shorter than that of the linear relaxation of the curved rolls. The total distance traveled 
by the boundary scales as e~^^^, where e is the reduced Rayleigh number. We obtain analytical 
I expressions for the relaxation rate of the modulation and for the time dependent traveling velocity 

^ . of the boundary, and especially their dependence on wavenumber. The results agree well with 

direct numerical solutions of the Swift-Hohenberg equation. We finally discuss the implications of 
I— I, our results on the coarsening rate of an ensemble of differently oriented domains in which grain 

f"]^ ■ boundary motion through curved rolls is the dominant coarsening mechanism. 
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I, . I. INTRODUCTION 

C ■ . 

Q ' This paper addresses the motion of a grain boundary separating two regions of parallel convective rolls, as can be 
O , observed in a Rayleigh-Benard convection cell of large aspect ratio. Each semi-infinite region (or grain) is comprised 
of roll patterns of same wavenumber qo, but the corresponding wavevectors are mutually perpendicular. Our main 
T-H , focus in this paper is the relationship between local roll curvature and grain boundary motion. Although the model 
^ ■ equations used in our analysis are appropriate for Rayleigh-Benard convection near onset Q, we expect that the 
[ qualitative features of our findings also hold in others systems that exhibit layered phases, such as, for example, 
■ lamellar phases in weakly segregated block copolymers ||J^. 

, Consider a disordered system brought into a layered phase, e.g. by a temperature quench in the case of a diblock 
' copolymer, or by a change in Rayleigh number in a Rayleigh-Benard convection cell. Domains comprised of rolls (in 
^-H ' Rayleigh-Benard convection) or lamellae (in block copolymers near a symmetric mixture) quickly form that have a 
• well defined characteristic wavenumber qo (in the case of Rayleigh-Benard convection near threshold, for example, 
qo lies on the marginal stability boundary against a zig-zag instability [^^). Due to translational and rotational 
invariance, the spontaneous relaxation from the initial disordered state leads in practice to a large number of such 
^ , domains, and a sufficiently large system remains isotropic macroscopically. Such a configuration also contains a large 

' density of defects, such as grain boundaries, disclination and dislocations. 
^ I Defect motion in two dimensional layered phases has been studied extensively, especially in connection with the 
O • evolution of convective rolls in Rayleigh-Benard cells . A primary question is how defect velocities are related to 
\ features of the background surrounding them, such as roll periodicity. It is well known, for example, that dislocations 
^ ■ climb 1^, or that grain boundary motion between two domains with straight rolls provides a wavelength selection 
k>( [ mechanism . On the other hand, less attention has been paid to the relationship between defect motion and roll 
; I curvature pi[ . We study in this paper a particular situation that appears to be prevalent (even if idealized) during 
the formation and evolution of layered structures: curvature induced motion of a grain boundary separating two semi- 
infinite ordered domains. We first show that if the rolls of one of the domains are periodically modulated along their 
transverse direction, the perturbation decays in time with a rate proportional to q*, where q is the wavenumber of the 
modulation. In addition, a novel feature is associated with this relaxation: the average position of the boundary does 
not remain stationary, but undergoes a net displacement such that the size of the region with straight unperturbed 
rolls increases. Analytic calculations based on the amplitude equation formalism show that this motion is a nonlinear 
effect, that it occurs in a time scale much shorter than the linear relaxation of the curved rolls. The total distance 
traveled by the boundary scales as e~^^^, where e is the reduced Rayleigh number. Consequently, the grain boundary 
can travel large distances, even for very small initial perturbations. 

The analysis of grain boundary motion and relaxation, while interesting in its own right, is also expected to 
contribute to our understanding of the coarsening of an ensemble of grains. Linear analyses of boundary or defect 
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motion have already been used in the past to predict coarsening exponents ||T2|,|3|. The method reproduces known 
exponents in the case of models A and B in the lexicon of Hohenberg and Halperin , and has also been used to 
predict coarsening exponents for 0(N) vector models jT^ (see also Coarsening laws for layered phases, on 

the other hand, remain largely unexplained ]r7| , p^ ,^|JT^-|20(| . 

In Section |l| we describe the configuration of the grain boundary studied, as well as the model equations. Section 



III presents a linear stability analysis of long wavelength transverse modulations near the boundary. The numerical 
method used to validate our solutions is outlined in Section IV, as well as the motivation for the nonlinear analysis 
given in Section Finally, in Section VI, we present conclusions, and briefly discuss possible implications of our 
results on coarsening rates. 



II. GRAIN BOUNDARY CONFIGURATION AND GOVERNING EQUATIONS 

The base or reference state is a stationary and planar grain boundary separating two semi-infinite domains A and 
B of straight convective rolls. Both sets of rolls have the same wavelength Ao — STr/go- It is in general expected from 
previous studies that such grain boundaries are stable against low wavcnumbcr perturbations |2l| ] . We wish to study 
here the decay of a perturbation of wavenumber g <C 90 applied in the direction transverse to the rolls of one of the 
domains, referred to as to domain A (see Fig. |l|a). In this article, we consider for all our analytic calculations the 
particular case in which the two sets of rolls are oriented at right angles relative to each other in the way depicted 
on Fig. IId. In this case, the rolls of domain A are parallel to the grain boundary itself. We expect that a 90° 
grain boundary is the boundary of lowest energy, and hence that it is the most common in an extended system that 
evolves spontaneously. However, our analysis can be generalized to cases such as shown on Fig. in which A and 
B rolls have arbitrary angles with respect to the grain boundary, and we do not anticipate qualitative changes in our 
conclusions. 

Our focus is on transverse modulations since this is the simplest kind of perturbation that induces roll curvature 
without any average roll compression (compression is associated with longitudinal perturbations along the roll normal) . 
Our study differs from that of ref. pTf in which a grain boundary was perturbed with a modulation of wave vector 
having components both transverse and longitudinal to the rolls. As was shown in ref. the longitudinal component 
of the perturbation dominated the relaxation. We analyze here the relaxation of a pure transverse (or curvature) 
mode, and obtain results that are qualitatively different from those of ref. |^ . 

Let an order parameter ijj represent, for example, the vertical velocity in the mid-plane of a convecting fluid layer. 
Just above onset, when the reduced Rayleigh number e = {R — Rc) / Rc 1 {Rc is the critical Rayleigh number for 
instability), ^ can be expanded as the superposition of two waves with slowly varying complex amplitudes A and B 

^P{x, y, t) = ^ [A{x, y, t) e^*^ + B{x, y, t) e'^"^ + c.c] ; A\,=^^ = 0, B\^^+^ = 0. (2.1) 

The amplitudes satisfy a set of two coupled Ginzburg-Landau equations that can be derived from the equation for tp 
(for instance the Swift- Hohenberg model |2|]) by multi scale analysis (|J^, 

dA^_6F dB^_SF 

dt SA ' dt SB' ^ ' ' 

where SF/SA is the functional derivative with respect to the complex conjugate of A, and F is a Lyapunov functional, 

F = Jdr^^ ldr{-e{\A\' + \B\') + ^{\A\^ + \Bf}+g^\A\'\B\' + 



- ^d'y)A\'+eo\idy - ,-dt)B\' . (2.3) 

The coherence length is of the order of q^^ ^ g and g±^ are interaction coefficients, and the time scale factor has been 
set to unity. These three parameters depend on the particular model equation considered for the original field ip. 



The nontrivial stationary solutions {Aq, Bq} of the governing set of equations (2.2)-(2.3) describe a planar boundary, 



and depend only on the coordinate x normal to the boundary. They are defined by, 

= eAo+eodlAo^gAl-g^B^oAo (2.4) 
= eBo- j^dtBo - gB^ - g±AlBo . (2.5) 
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This system of equations was extensively studied in refs. ||9|,y^. It is a planar grain boundary of thickness proportional 
to Co/'s/s- Note that the system of equations is not invariant under permutation of Aq and Bq. The amplitude Aq 
of the rolls parallel to the interface vanishes as exp{xy/e/^o) when x —oo, and saturates to (e/g)^/^tanh(x-ye/^o) 
when X —^ +00. The behavior of the amplitude of the rolls perpendicular to the interface is slightly different: 
Bq{x) — (e/g)^/^ cx exp(a;Y^/^o) when x —00, and there is a point x* such that Bo{x*) — 0. It is customary to 
define the location of the grain boundary by the point x*. To a good approximation, Bq{x) ~ for x> x* . 

We emphasize that a stationary solution only exists when the wavenumber of the solution equals the marginal 
wavenumber go, and that at first order in e, the solutions Aq and Bq only depend on the slow spatial variable e^/^x, 
but not on the fast scale x. Hence, the location of the boundary and the phase of the rolls are independent H. 



III. LINEAR STABILITY ANALYSIS OF A 90° GRAIN BOUNDARY 



We linearize Eqs. (2.2) around {Aq, Bq}, and assume the following perturbed solutions (note that the form of the 
perturbation explicitly assumes that the boundary does not undergo a net average displacement), 



A{x,v,t) = {AQ{x) + ~a)e'^, 
B{x,y,t)=Boix)+b . 

The phase (/) represents a transverse distortion, as shown in Fig. |l]b, and is of the general form 



(/) = 0(a;) cos(gj/)e° 



< 1. 



The real fields a and b are amplitude corrections and similarly read, 

a — a{x) cos(gy)e'^* 
b = b{x) cos{qy)e'^* . 



(3.1) 
(3.2) 



(3.3) 



(3.4) 
(3.5) 



In the above expansion, we have neglected the imaginary part of B. Im{B) can written as dix) sm{qy)e'^* ; yet, 
because of the fourth order x-derivative in the equation satisfied by B (the second relation in Eq. ( |2.2|)), d{x) is of 
order e lower than the real part function b{x), and hence negligible. Inserting Eqs. (3.1)-(3.2) into Eqs. (|2.4), we find. 



aa = ea + ^lidla - {qV^ll)a ~ (gV9o)5.(<^^o)] - {igAl + g^Bl)a - ig^A^B^b 
ab = eb- ilq^b - [igBl + g^Al)b - ig^B^A^a . 



(3.6) 
(3.7) 
(3.8) 



We further expand the solutions in power series of the small parameter , and define, a — ao + q^a2 + q'^a-A + •• 



-q^(t>2- 



q 4>i + .. ; 6 = feo + q b2 + q bi + ... Since the limit q 



corresponds to a uniform translation of the 



interface, a = q <72 + f? cr4 + •.- We now analyze the resulting equations order by order in q . 



At order 5°, Eq. (3.7) reduces to Aldx<po = est. Since Aq{x) ^ diverges exponentially at —00 and dx4> must remain 



finite on the whole interval [—00,00], the only admissible solution is dx4>Q = 0. Hence (j)o =cst is a free parameter 
that represents the magnitude of the initial phase modulation. With the notation of Fig. |i|d, (f)Q is simply related to 
the magnitude of the roll deformation 6x0 through 0o = qoSxo {6x0 <C Aq). At order Eqs.(|3.6|) and (3^) can be 
written as 



Hn 



ao 
60 



= 0. 



If we take the a;-derivative of Eqs. (^^) and (|]^) we recover the above equation, with oq and bo replaced by dxAo 
and dxBo- Hence the solutions at this order are: 



ao{x) = aodxAo, bo{x) = aodxBo, 

where ao is a constant to be determined from the solvability condition at the next order. 
At order q^ , the solutions 02 and &2 satisfy, 



Ho 



02 
62 



(aoCT2 + 4io£.l/ qn)dxAo 
ao(CT2 + £.o)dxBo 



(3.9) 



(3.10) 
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The solvability condition at order requires that the right hand side vector be orthogonal to the kernel of the adjoint 
of iJo- Since Hq is hermitian, and HoidxAo, d^Bo) = 0, we find, 



(aocra + Mo/lo) / (d^Aofdx + ao{(T2 + Co) / {dxBofdx = 



(3.11) 



On the other hand, at order one can replace a4> by a24>o, and a by a^dxAQ in Eq. ( |3.7|) . After multiplying Eq. 
(3.7) by Aq, and integrating over x, one finds 



Aldx(t)2 = 



so 



A^[u)du — 



90 



Ao(x)a^Ao(a;) - / {duAaYdu 



(3.12) 



where C is a constant of integration. However, the condition that |92;02| < oo when x — > —oo requires C — Q. 
Furthermore, in the limit x — > +cx), ~^ (^/s)^^^; hence A^{u)du oc x. The gradient of 02 remains finite only if 
(72 = 0. Therefore Eq. ( |3.11 ) with (T2 = now yields the value of the solvability constant ao, 



-5xq 



(dxAofdx 



{dxBofdx 



(3.13) 



where the definition (pQ = qoSxo has also been used. 

To summarize our results up to this point, the perturbed amplitudes at order g*^ are given by, 

A = [Ao{x) + aoCOs{qy)e^*dxAo{x)] e'^ ~ Ao (x + ao cos(gy)e'^*) e'ioS.ocosigy)e'* ^ 
B Bq (x + ao cos((7y)e'^*) , 



(3.14) 
(3.15) 



where we have used Eqs. (p.l|)-( |3.2| ) and ( |3.9| ). Equations ( 3.14 ) and ( 3.K ) show that the amplitude moduli \A\ and 
of the weakly modulated rolls at this order equal the unmodulated profiles with a simple change in the coordinate 
origin. By analogy with the planar c ase |l9| , |lO| , we define the location of the modulated boundary by the set of points 
{xg, y} such that B{xg, y) = 0. Eq. (3.15) indicates that the location of the boundary is given by the curve 



Xg(y) = X* - ao cos(gy)e'^* , 

with i?o(x*) = 0. It is a local translation relative to the planar boundary by — ao cos(q?/)e'^* along the x-axis. Note, 
on the other hand, that the lines of constant phase of A are given by x = est — 5xq cos(gy)e'^* instead, so that 
— (5xo cos((72/)e'^* represents the local deformation of the straight rolls. As shown by Eq. ( |3.13| ), boundary and roll 
deformations are not independent: the lengths ao and 5xq are related to each other through the non- uniform profiles 
Aq and Bq. A qui te non-intuitive result is the opposite directions of the related translations, evidenced by the minus 
sign in Eq. ( 3.13 ). Consider for instance a point where the lines of constant phase of A are displaced towards the 
region B] the actual position of the boundary, however, is displaced towards region A. This effect can be seen more 
easily on Figure || (see also Section IV). 

At order q'^ integration of Eq. (3/7) leads to. 



Ai^dx(i 



1± 



1 



(j)o^{u)du - 



1 



90 J- 



Ao{u)dua2du, 



where the second order result (T2 — has been used. The second integral in the right hand side is finite, but the first 
one diverges as x when x — > +oo, except if 0-4/^0 + l/4<Zo ~ 0. We therefore conclude that, at leading order in q 



'Ml 



(3.16) 



This is one of the central results of this section: the modulated boundary is expected to relax exponentially with a 
rate proportional to g^. 

The linear analysis shows that the interface is always stable with respect to long wavelength perturbations. This is in 
agreement with the case studied in ref. j2^, although the authors concluded that the decay rate a (x q^ instead. Such 
a behavior was found because the wave vector of the modulation considered in their work had a nonzero component 
in the longitudinal direction. We note that the decay rate of Eq. (3.16) is in fact identical to that of a single wave 
(A-roUs only in Figure [l]) 0. According to Eq. (3.16), the relaxation of the boundary is completely determined by 
the relaxation of the A rolls far from it. 
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Unfortunately, in spite of the fact that the main predictions of the perturbation analysis agree well with a numerical 
solution of the Swift-Hohenberg equation, the linear analysis just discussed is not uniformly valid (we will further 
elaborate on this point in Section^. Consider Eq. (3.12) for <^2 in the limit x — > +00. The leading behavior is given 
by 



qaAlioo) 



{duAofdu X . 



(3.17) 



Therefore the phase perturbation 



- q'^4>2, which is assumed to be small in the derivation of Eq. (pj), diverges 



as X for large x through the second order correction 
y-component of the wavevector of the roll pattern is. 



at any finite time t. As a consequence, from Eq. (|3.3D, the 



dy4> ~qsiii{qy) e'"*((/)o + g^(/)2), 

and also diverges with x. This implies that, far enough from the boundary, the y-component of the wave vector 
becomes larger than the x component, given hy qo -\- dx(t) — qo- This solution represents an increasingly zig-zagging 
wave, which is of course not physical. As it is apparent from Eq. (3.17), the singularity arises because of the presence 
of a non- uniform amplitude of the pattern near the boundary O^Aq ^ 0. In fact, we think it is possible that the 
linearized system with non-constant coefficients (3.6)-(3.8) does not have bounded eigenfunctions. 



It appears that the failure of the linear analysis is related to an additional important feature of boundary motion, not 
taken into account in the previous calculation: a net translation of the average position of the boundary towards the 
A-region, as shown by the numerical calculations presented in the following Section. After we present the numerical 
evidence, we argue in Section ^ that this feature cannot be obtained at the linear level in the amplitude of the 
perturbation, and that it is in fact a singular contribution in the limit e ^ 0. 



IV. NUMERICAL SOLUTION OF THE SWIFT-HOHENBERG EQUATION 

In order to test the predictions of the previous Section, and to obtain further insight into grain boundary motion, 
we have numerically solved the Swift-Hohenberg model of Rayleigh-Benard convection, 



5^ = ^^ - 4^ 



1^ = eV' - 73 (9o + V^)> - ^P'. (4.1) 



We consider the evolution from an initial condition that corresponds to the geometry shown in Fig. |]b. We note 
that solving Eq. ( [4.l| ) in the case of curved rolls is much simpler than solving the two dimensional Ginzburg-Landau 
equations ( |2.2|). 

Equation ( [l.l| ) is solved numerically with a pseudo-spectral method. Further details on the algorithm and the time 
integration scheme can be found in ref. . The stability of the algorithm allows relatively large values of the time 
step, which is fixed to 0.4 in the dimcnsionlcss time units of Eq. (4.1). Equation ( [4.1| ) is then discretized on a square 



grid of size 512 x 512, and occasionally 256 x 256. Spatial discretization is such that there are 8 grid nodes per 
wavelength Aq. In all the following numerical examples, we have chosen ^0 = "^Qo^^ which is close to the value that 
corresponds to the stress free boundary conditions in the original system of fluid mechanical equations, from which the 
Swift-Hohenberg equation is derived as an approximation. We use periodic boundary conditions in both directions, 
and hence the initial condition comprises two symmetric, well separated grain boundaries located at a; = 1/4 and 
x = 3/4, in units of the system size. We have verified that the numerical solutions of the Swift-Hohenberg equation 
for the stationary straight grain boundary coincides with those obtained by directly solving the amplitude equations 



(2^)- (2^), with the appropriate parameters that follow from the multi-scale analysis of the Swift-Hohenberg equation 
(5 = 3/4 and g± — 3/2). The modulated initial condition is implemented with the help of the one dimensional solution, 
as 

tfjix, y,t = Q)= Ao{x) cos[go2; + qoSxQ cos{qy)] + Bo{x) cos[qQy]. (4.2) 

A typical configuration is shown in Figure H, taken at the intermediate time t = 1600, with e — 0.04, q = qo/l6 
and a system size of 256 x 256. Note that this solution is consistent with our earlier assumption of neglecting Ira{B) 



in the perturbation equations (Eq. (3.2)), since the B rolls remain straight near the grain boundary despite of the 



transverse modulation of the A rolls. The opposition in phase of the roll profile and of the boundary location, given 
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in Eq. (3.13), is also well reproduced numerically Yet, this effect is easy to observe only for relatively large values of 
Sxq, as shown on Fig. ^j. 

Figure ^ shows our numerical results for the inverse decay rate —a~^ as a function of q. These results are obtained 
for a system of size 512 x 512, so that the number of rolls that separate the two grain boundaries is twice larger 
than in Figure |2| We checked that the planar grain boundaries did not appreciably interact, and therefore remained 
stationary. The value of a was determined from an exponential fit to the decay of the phase of a given roll, for a 
small initial Sxq- We first compute numerically a in the absence of any grain boundary, i.e. with the single wave A. 
We find very good agreement with Eq. ( ^.16 ), given that there are no adjustable parameters in the theoretical curve. 
When grain boundaries are present, the numerical results (square symbols) are somewhat higher, but still compare 
well with Eq. (3.16). They remain closer to the law a oc —q"^, rather than to, say, a oc —q^, and certainly than to 
a oc —q^. The small discrepancy observed is probably due to the finite size of the system: although the two grain 
boundaries are separated by roughly 30 rolls, they slightly interact during relaxation. As a check, we have computed 
a{q = 3go/32) again in a system of size 1024 x 1024, with same value of Aq. The difference between that numerical 
result and the theoretical curve is reduced by a factor of 2.5 compared with a system size of 512 x 512. 

The numerical solution reveals an additional feature which was not included in our linear analysis: while the 
boundary modulation is relaxing, there is net motion of the grain boundary towards the A roll region. This motion is 
a pure amplitude mode, as the lines of constant roll phase remain stationary on average. Therefore, after a long time, 
the modulated rolls have completely relaxed but the size of the region of A rolls has decreased, and the grain boundary 
is located at a distance d^o away from its initial average location. The distance dao exceeds the initial deformation 
Sxq, and it can be even significantly larger than the roll wavelength Aq. In the configuration shown in Fig. || (system 
of size 256 x 256), the magnitude of the perturbation was deliberately chosen relatively large {6xo = Aq), so that, at 
late times, the modulated region completely disappears by annihilation of the two grain boundaries. 

We now report our results for doo in a bigger system of size 512 x 512, and with q = 3go/32. Figure ^ shows our 
results for the total traveled distance doo as a function of 5xq. The complete relaxation of the perturbation requires 
times around 10^, and the location of the (nearly) flat boundary is again defined to be the point x* at which B{x*) — 
This point if obtained from the field -0 through the relation B{x) = [i^{x, y = Aq) — 'ipix, y = 3Ao/2)]/2 for a flat 
interface (see Eq. (2.1)). We observe in the figure that the total traveled distance increases nonlinearly with 5xq. The 
following Section gives an interpretation of this feature, and further discussion of the results. 



V. DEFECT MOTION THROUGH ENERGY RELAXATION 



We study in this section a more general form for the perturbed amplitudes than those given in Eqs. (3.14)-(3.15) 
We assume a perturbed solution of the form, 



where 



A^Ao{x- K5x{t) cos{qy) - d{t)) e**«^(*)™''(9y) 
B = Bq {x — KSx{t) cos{qy) — d{t)) 



5x{t) = 5xoe'"*, a = --^ (7^ and k = / {d^Aofdx/ / {d^Bofdx 



(5.1) 
(5.2) 



(5.3) 



The quantity d {t) r epresents the net distance traveled by the grain boundary at time t. A linear stability analysis 
based on (pj|)-(5.2) along the lines of Section [D leads to d{t) = for an interface initially located at the origin. As 
we show below, the net displacement depends nonlinearly on the initial modulation, and hence cannot be obtained 
from a linearized set of equations. We use here a different method that has been widely used to study of defect motion 
in potential systems jl) . It is based on the identity. 



dt 



dr 



dA 
'dt 



dB 
dt 



(5.4) 



directly derived from Eqs. (2.2)-(2.3) 



The left hand side of Eq. (5.4) can be calculated approximately by substituting Eqs. (5.1) and (5.2) into Eq. (p.g|), 



and taking the time derivative. After some algebra, we find. 



dF_ 
lit 



— = I dr\ -d dxTo 



SO? 



cos^ (qy) 26x5xAI{u) — Sx^d dxAl{u) 



(5.5) 
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up to order 5x^ . is the free energy density of a planar boundary, and u = x — K5x{t) cos{qy) — d{t). The term 
involving J^q vanishes since J^o{x = +00) = To{x = —00). In Eq. ( |5.5D , we have neglected the contributions from 
the terms dyAo^u) and dyBo{u); they are of order e^/^ smaller than Ao^ujdye^'^ because of the slow variations of the 
amplitudes compared with the roll periodicity. In addition, we have used the approximation dtAo{u) ~ —d dxAQ{u) 
by neglecting —KSx{t) cos{qy) dxAQ{u): the contributions proportional to this last term vanish with J cos{qy)dy at 
leading order in e when spatial integration over y is performed in Eq. (5.5). For the same reason, all terms proportional 
to Sx in Eq. (5.5) do not contribute. 

By using similar considerations, the right hand side of Eq. (p^) at leading order in e and Sx{t) reads, 



df{d^[idMu)f + {d.Boiu)f] + idx)^q^o cos2(gy)A2(u 



'} 



(5.6) 



By combining Eqs. { p.5\) and (|5.6|), and by using Eq. ( |5.3| ) for Sx{t), we find that all terms that do not involve d 
cancel, and obtain 



m = ( 4?' ) 77 



[qoSx{t)]^ 



Mr J 4g j^^dx[id,Aor + id,Bor] ' 



(5.7) 



where we have used the identities Aq{+oo) — e/g and dy cos^ (qy) / L = 1/2, where L is the length of the grain 



boundary. Equation (5/7) shows that d > 0; i.e., in the configuration shown in Fig. |l| the motion of the boundary is 
such that the modulated rolls are progressively replaced by the advancing perpendicular straight rolls. Furthermore, 
and contrary to the usual motion of kinks, the velocity is not constant but decreases as Sx{t)^. Since at long times 
Sx{t) ~ 0, the boundary eventually stops. In Appendix A, we show that Eq. (5/7) generalizes the equation describing 
the uniform motion of a grain boundary through a time independent modu late d background. The method presented 
in Appendix A also provides a simpler alternative way to extend the result (5.7) to the case in which the transversally 
modulated rolls are not parallel to the grain boundary. 

Equations (5.3) and (5.7) can now be used to calculate the total distance traveled by the grain boundary doo, 
allowing further comparison with the numerical results of Section LV. We find, 



doc — ^ 



{qaSxoY^ 



8ff JZodx[{d.Aor + idM^^ 



(5.8) 



Note that, in contrast to d, d^o does not depend on the wavenumber of the modulation q, for small q. Since 
Aq{x) — yJTJgf{^/e x/^q) at leading order in e, and if a similar form is assumed for _Bo, it is easy to see from Eqs. 
( ^.7|) and (5^) that d and doo are proportional to e^^/^ |^^. Hence th e dis place ment of the boundary diverges close 
to onset, regardless of the initial amplitude of the modulation. Eqs. (5/7) and ( |5.8|) also show that the magnitude 



of this displacement depends quadratically on the perturbation 5xq, explaining why it can not be derived from the 
linear analysis presented in Section III. Finally, we note that the velocity d{t) is proportional to t~^^'^5x{t)'^q'^^ at 
leading order in e, 5x{t) and q. 

The mechanism for net grain boundary motion can be understood as the follows: Relaxation requires energy decrease 
while straight rolls have the same energy regardless of their orientation (both base states A and B, are in principle 
degenerate). However, when the boundary is modulated in the way depicted in Fig. ^ only the rolls parallel to the 
boundary are distorted at lowest order in e, whereas the perpendicular rolls remain straight. Therefore, the energy 
of the modulated configuration can decrease by both linear relaxation of the initially curved rolls and by nonlinear 
net displacement of the grain boundary, the effect of which is to replace curved rolls by straight ones. Close to onset, 
relaxation through the nonlinear mode is dominant, and even a small perturbation 5xq can induce a large boundary 
displacement doc ^ 6xq. Of course, for fixed e, doo goes to zero with 5xq. It is precisely the observation that the 
amplitude of this nonlinear mode diverges at e ^ that suggests that the linearized problem discussed in Section III 
may not have any bounded eigenfunctions. 

At fixed e, it is instructive to obtain the magnitude of 5xo above which the velocity of the boundary is much larger 
than the phase velocity of the rolls; i.e., 

rf(t = 0) > |5i(t = 0)|. (5.9) 

In this limit, boundary motion occurs at approximately constant roll phase. Given that 5x{Q) = adxQ, and by using 
Eq. (|^ for d(0), Eq. ^) leads to 

2 



Sxo » (^^^ ^ J°° dx [{d^Aof + id.Bof] 
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With the numerical solutions obtained for Ao{x) and Bq{x) from the coupled Ginzburg-Landau equations, we find 
that Sxq 3> 5.0 10~^Ao for e = 0.04, and 5xq ^ 7.4 10~^Ao for e = 0.1. Therefore we conclude that nonlinear motion 
dominates linear relaxation even for relatively small initial perturbations. We recall that our perturbative calculation 
is valid if 5xq < Aq. 

We have verified that dao does not appreciably depend on the wavelength of the modulation, as shown in Eq. 
(5.8). Figure ^ shows our analytical and numerical results of the total distance traveled by the grain boundary as a 



function of 5xq. The theoretical curve has no adjustable parameters, and is based on the numerical resolution of Eqs. 



(2.4)-(2.5) with the values g = 3/4 and g± — 3/2 that correspond to the Swift-Hohenberg equation. We performed 
two series of runs, for e = 0.04 and e = 0.1. At e = 0.04, Eq. ( |5.8| ) is in quantitative agreement with the numerical 
data for small modulations. The discrepancy observed at large modulation indicates the breakdown of the linear 



relaxation ol 5x{t): perturbations start decaying at a faster rate than (B.16), so that Eq. (5.8) is over-estimating the 
traveled distance at short times. At e = 0.1, we expect a smaller displacement doo, for fixed 5xq, than in the case 
e = 0.04. This is what we observe. However, the numerical results show that d^o does not vary continuously with Jccq, 
but rather it increases discontinuously. As it is apparent in Figure ^ (crosses) , doo takes now only multiple values of 
half the wavelength of the pattern. A lower threshold is needed to displace the grain boundary, which we estimate to 
be 5xc — Ao/7. These observations cannot be explained by the present analysis, and reveal non-adiabatic effects that 
go beyond the amplitude equation formalism. Such effects have been previously observed in other studies of fronts 
propagation [p]j25|]. The underlying cause is that with increasing e the variations of A in the slow variable e^/^x do 
not decouple with the fast variable describing the underlying periodic structure. In this case, the solution for the 
amplitudes depends on the position of the rolls. 

VI. CONCLUSIONS 

We have found an additional mechanism for grain boundary motion in layered phases which is driven by roll 
curvature in the vicinity of the boundary. We have analytically studied the motion of an isolated grain boundary 
within the amplitude equations formalism, and numerically by direct solution of the Swift-Hohenberg equation. 
Transverse perturbations of the rolls near the grain boundary decay exponentially in time with a rate proportional 
to q^, where q is the wavenumbcr of the perturbation. In addition, we have found a nonlinear mode involving net 
translation of the average location of the boundary. On the basis of functional minimization, this mode is preferred 
to simple linear relaxation near onset. The net velocity of the grain boundary has been computed, and is given 



in Eq. (5.7). This result can be interpreted as giving the velocity as the ratio between a time dependent external 
force imposed by curved rolls, and a drag term or mobility that depends on the local amplitude profile of the grain 
boundary. Close to onset, the drag term goes to zero as e^/^, and makes the interface velocity (and the total traveled 
distance) diverge as e~^/^. We have argued that this motion cannot be explained by linear analysis. 

The precise relationship between these results and the related problem of coarsening of layered phases requires 
further study. However, we can offer a few comments based on the results presented in this paper. If the temporal 
evolution of a disordered configuration were controlled by the relaxation of an ensemble of grain boundaries moving 



through a background of curved rolls, then Eq. (5/7) could be used to infer the coarsening rate. If w is a c har acteristic 



grain boundary speed, and k, is the characteristic curvature of the rolls ahead of the boundary, then Eq. (5.7) leads to 
V (X . If self-similarity holds during the coarsening of the structure, this last relationship would imply a coarsening 
law l{t) ^ t^/^ , where l{t) is any measure of the linear scale of the structure. This result for the coarsening law 
disagrees with previous literature on the subject [p7| , p"2P ,[r^-po| . 

Although the configuration studied in this paper is idealized, is is conceivable that it describes one among possi- 
bly many competing mechanisms during coarsening. A numerical solution of the two dimensional Swift-Hohenberg 
equation close to onset shows that curved rolls are essentially immobile due to topological constraints (mostly discli- 
nations), whereas grain boundaries move over large distances. It is likely that the motion of the latter is in part due 
to a background of curved rolls with a characteristic curvature which is set by the spatial distribution of disclinations. 
Disclinations, in turn, can be eliminated by grain boundary motion. If the characteristic length scales associated 
with both defects (grain boundary perimeter and disclination separation) are proportional to each other, as required 
by self-similarity, then t^/^ would be a possible contribution to the coarsening law. Whether this mechanism would 
dominate the asymptotic behavior as I ^ oo \s being currently investigated. 



8 



ACKNOWLEDGMENTS 



This research has been supported by the U.S. Department of Energy, contract No. DE-FG05-95ER14566, and also 
in part by the Supercomputer Computations Research Institute, which is partially funded by the U.S. Department of 
Energy, contract No. DE-FC05-85ER25000. 



APPENDIX A: GRAIN BOUNDARY MOTION WITHOUT PHASE RELAXATION 



Equation (5.7) can be derived in an different way, which is similar to the calculation of the velocity of a climbing 
dislocation [y! Suppose that a set of rolls is modulated along the transverse direction, with wave vector q, and is 
rigidly held with a time independent modulation amplitude 5xf. Hence, the system is allowed to evolve by rigid 
translation alone, and no roll (or boundary) relaxation is allowed. We then introduce the perturbations (5.1)-(pT^) 
with 6x{t) — 6xf. Equation ( ^.4| ) can be simply recast as 



d J dy{T{x - -oo, y) - T{x = +00, y)} = -2d^ J dr [\d.,A\^ + \d,B\^] 



(Al) 



where d is the constant velocity of the grain boundary. The left-hand side of Eq. (Al) represents the free energy lost 
per unit time by the system as the grain boundary moves. The straight rolls, of lower free energy, replace curved 
ones. The difference dF is computed that corresponds to increasing the area of the domain occupied by B straight 
rolls by an amount d dt, and by decreasing by the same amount the area of the domain occupied by A curved rolls. 
The free energy gain is easily computed by using Eq. ( p.3| ) with B{—oo) = {e/g)^^^, A{—oo) = 0, and B{+oo) = 0, 
A{+oo) = {e/g)^/'^e''"> ^^f^°^iiv). In the right-hand side of Eq. (^), \dxA{B)\ is replaced by \d^Ao{Bo)\ at leading 
order. Eq. ( |AlD yields Eq. (5.7) straightforwardly, where 6x{t) must be replaced by Sxf. By using Eq. ( |AlD , the 
numerator of Eq. (5.7) can now be interpreted as a (time dependent) external force acting on the line defect, the 



denominator being the drag term or friction coefficient. It is interesting to note that the more rigorous analysis of the 
combined effects of roll relaxation and front advance that leads to Eq. (^^) gives the same expression as when roll 
relaxation is omitted. Beside being simpler, an additional advantage of the energy argument presented here is that 
it can be applied to a configuration where the set of rolls A and B are not parallel nor perpendicular to the grain 
boundary. In this case, the external force J dy{J-'{x = —oo,y) — T{x — +00, y)} remains unchanged, by rotational 
invariance. Only the drag term would probably differ, because of the different grain boundary profile. However, we 
expect that the scaling behavior d oc 5x^q^ still holds. 
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FIG. 1. a) Grain boundary separating two set of rolls A and B of same periodicity {\qo\ — \%'\ = qo). The rolls of domain 
A are weakly curved by a transverse modulation of wavenumber q ^ qo- b) Particular case studied here, corresponding to a 
90° orientation. Sxo is the magnitude of the phase modulation. 



FIG. 2. Configuration with curved interfaces obtained by numerical solution of the Swift-Hohenberg equation in a square 
grid with 256 x 256 nodes, with e = 0.04, q = qo/16 and Sxo = Ao = 8. Because of the curved rolls, the grain boundaries are 
moving towards each other. 



FIG. 3. Decay time —a ^ (in dimensionless time units) of the grain boundary modulation as a function of wavenumber 



q. The square symbols are the numerical results while the solid line corresponds to Eq. (3.16). 
numerical results obtained from the modulation of a single wave, without any grain boundary. 



The cross symbols are the 
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FIG. 4. Total distance traveled by the interface dao as a function of the magnitude of the initial modulation Sxq, for e — 0.04 
and e — 0.1. {q = 3go/32, the grid size is 512 x 512 and Ao = 8). The symbols are the numerical results, and the solid line 
corresponds to Eq. (5.1 
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